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ABSTRACT 

Ultra-high-energy cosmic rays (UHECRs) may be produced in active galactic nuclei 
(AGN) or gamma-ray burst (GRB) jets. I argue that magnetic reconnection in jets can accel- 
erate UHECRs rather independently of physical processes in the magnetic dissipation region. 
First order Fermi acceleration can efficiently take place in the region where the unrecon- 
nected (upstream) magnetized fluid converges into the reconnection layer. I find that protons 
can reach energies up to £ ~ 10^^ eV in GRB and powerful AGN jets while iron nuclei can 
reach similar energies in AGN jets of more moderate luminosity. 

Key words: Cosmic rays: ultra-high energies - gamma rays: bursts - galaxies: active - mag- 
netic reconnection 



1 INTRODUCTION 

The origin of UHECRs (particles with energy £^^10^^ eV) remains 
a mystery. The fact that magnetic fields cannot confine them in the 
Galaxy and the observed spectral cutoff at ~ 6x10^^ eV point to ex- 
tragalactic sources at these energies (Abbasi et al. 2008; Abraham 
et al. 2008). For any astrophysical source to accelerate particles to 
the highest observed energies of Em^^O^^ eV, tight constraints need 
to be satisfied. Possibly, the tightest come from the need for parti- 
cles with energy Em to be confined within the size of the source 
(Hillas 1984) and the acceleration to take place within the available 
(dynamical) time of the system. Both constraints may be met in 
very powerful sources. Viable astrophysical sources for UHECRs 
are relativistic jets from active galactic nuclei (see,. e.g. Halzen & 
Hooper 2002) and gamma-ray bursts (Milgrom & Usov 1995; Wax- 
man 1995;Vietri 1995; see, however, Ghisellini et al. 2008; Inoue 
2008 for alternatives). 

The acceleration mechanism for UHECRs is usually postu- 
lated to be the first order Fermi mechanism in (mildly) relativistic 
shocks in the jets (because of internal or external interactions; Gal- 
lant & Achterberg 1999; Achterberg et al. 2001). For the UHECRs 
to be confined by shocks both the upstream and downstream need 
to be strongly magnetized (i.e. with magnetic energy density not 
much less than the total energy density; e.g. Waxman 1995). The 
magnetic fields cannot be of small scale (generated by plasma insta- 
bilities) since in this case the small angle scattering of the energetic 
particles makes the acceleration process too slow to be of relevance 
(Kirk & Reville 2010). Strong, large scale fields are also unlikely 
to work since they do not allow for the particles to repeatedly cross 
the shock front (where the acceleration takes place) for most field 
inclinations (e.g. Sironi & Spitkovsky 2009). 
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Here I discuss an alternative mechanism for UHECR accel- 
eration in a relativistic flow. This mechanism applies to magnetic 
reconnection regions in Poynting-flux dominated flows (e.g. flows 
with Poynting-to-kinetic flux ratio cr > 1). The particle accelera- 
tion takes place through the first order Fermi mechanism because 
of particles reflected in the magnetized plasma that is converging 
in the reconnection region (see Speiser 1965 and de Gouveia dal 
Pino & Lazarian 2005 for a similar mechanism applied to Earth's 
magnetotail and the ejections from the microquasar GRS 1915+105 
respectively). The strong, large scale field in the upstream of the re- 
connection region that approaches the reconnection layer at subrel- 
ativistic speed results in a very efficient acceleration configuration. 

In the next section, I describe the mechanism for particle ac- 
celeration to ultra-high energies. This mechanism is applied to rel- 
ativistic jets in Sect. 3. Discussion is given in Sect. 4. 



2 PARTICLE ACCELERATION IN MAGNETIC 
RECONNECTION REGIONS 

In this work, I do not assume any specific mechanism for mag- 
netic reconnection. I consider the rather generic geometry where 
magnetized fluid with reversing polarity over a lengthscale /rec is 
advected into the dissipation layer with speed /?rec = Vrec/c- The 
plasma is heated/compressed and accelerated by the released mag- 
netic energy and leaves the region through a narrow layer of thick- 
ness 6 <^ /rec at the Alfven speed of the upstream plasma /?out ~ /^a 
(e.g. Lyubarsky 2005; see fig. 1). For high <t flows considered here 

/5out - 1. 

Particle acceleration at the reconnection layer can be compli- 
cated by the exact, small-scale reconnection geometry. At the loca- 
tion where magnetic energy is dissipated one may deal with current 
sheets (Parker 1957; Sweet 1958), slow MHD shocks (Petschek 
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1964; Zenitani, Hesse & Klimas 2009), MHD turbulence (Lazarian 
& Visniak 1999; Matthaeus & Lamkin 1986; Loureiro et al. 2009) 
and/or secondary tearing instabilities (Drake et al. 2006; Loureiro, 
Schekochihin & Cowley 2007; Samtaney et al. 2009). Regardless 
of the details of the reconnection mechanism, however, any suc- 
cessful reconnection model must explain the observed fast rate with 
which magnetic field lines are advected into the reconnection re- 
gion (e.g. Lin et al. 2003). The so-called reconnection speed has to 
be a substantial fraction (a tenth or so) of the Alfven speed p^ of 
the upstream plasma. For^^A ~ 1 the reconnection speed is subrel- 
ativistic: /5rec = 6 ~ 0.1. 

I assume that a preacceleration mechanism is able to acceler- 
ate ions to large enough energy so that their gyroradius Rg becomes 
larger than the thickness of the layer 6. Such preacceleration may 
take place in contracting magnetic islands (Drake et al. 2006), cur- 
rent sheets (e.g. Kirk 2004) or slow MHD shocks. After the initial 
acceleration (discussed in Sect. 3.3), the particle trajectories and 
energy will be subject to at most mild change while they cross the 
dissipation layer. Still, as long as Rg < /rec, the particles are confined 
in the wider reconnection region and can be accelerated further. 



2.1 The acceleration cycle 

The acceleration of the particles in the reconnection region may 
be viewed from two equivalent perspectives. In the rest fame of 
the current layer, the upstream region contains an advective elec- 
tric field directed along the x axis of strength fix = PxqcB (the 
so-called reconnection electric field). A particle bounces back and 
forth around the reconnection layer in a betatron-like orbit (also 
called Speiser orbit) schematically shown in Fig. 1. It is continu- 
ously accelerated by the electric field with its energy increasing lin- 
early with distance x traveled along the layer E' ~ ep^ecBx (Speiser 
1965). The acceleration may also be viewed as result of repeated 
magnetic reflections in the upstream flow. The fact that the up- 
stream is converging towards the dissipation region, allows for a 
first order Fermi acceleration to operate (de Gouveia dal Pino & 
Lazarian 2005). As a result, particles gain a fixed fractional energy 
per cycle that can be estimated by the following considerations (that 
involve two Lorentz boosts). 

Consider that a particle of Lorentz factor 71 » 1 (measured 
in, e.g., the lower side of the upstream region; see Fig. 1) enters 
at an angle 6 (with respect to the normal to the reconnection layer, 
i.e., the y axis) into the opposite side of the upstream. The relative 
speed of the two regions is 



/^r- 



1+/ 



(1) 



The Lorentz factor 72 of the same particle measured in the rest 
frame of the upper side is given by the transformation 72 = 7i7r(l + 
Pj: cos 6). The particle performs a fraction of a gyration and leaves 
the upper side at an angle -^[jThe particle completes the cycle by 
returning to the lower region with Lorentz factor (in the rest frame 
of that region) 73 = 727r(l + Pr cos 6). The amplification in energy 
during the 1 ^ 2 ^ 1 cycle is 



^ In general the particle will return at an angle l^'l different than 6. The 
change in the angle is small during a single cycle for subrelativistic /?rec 
discussed here and can be ignored. Furthermore, a small deflection of the 
particle may take place while it crosses the reconnection layer that intro- 
duces some degree stochasticity to the orbit but does not affect the main 
arguments presented here. 
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Figure 1. Sketch of the reconnection geometry. High-cr plasma approaches 
the reconnection layer at subrelativistic speed ySrec from both sides (up- 
stream flow). Plasma leaves the region through a thin layer with ySout ~ 1 
perpendicular to the plane of the sketch (downstream flow). Energetic ions 
are Fermi accelerated by been repeatedly reflected by the upper and lower 
sides of the upstream. 
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In the limit of ^ = (particle crossing perpendicular to the 
reconnection layer), the amplification A is given as function of /?rec- 



A(0): 
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(3) 



where eq. (1) is used in the last step. If any scattering process 
(e.g. when crossing the reconnection layer) keeps the particle dis- 
tribution quasi-isotropic, particles will cross the dissipation region 
with a angular probability distribution P(6)d6 = 2 sin ^ cos 6d6 with 
0^6^ 7t/2. In this case, the average amplification per circle is 



<A> 



nn/2 
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A(e)P(6)d6 = y^[\ + ^P, + ^P^l 



1 



(4) 
73 as found in de 



Note that, in the limit /?rec ^ 1, < A >= 
Gouveia dal Pino & Lazarian (2005). 

In Fig. 2, 1 plot the A(0) and < A > as a function of /?rec- Note 
that A(0) and < A > have rather similar values making the results 
that follow rather insensitive to the exact angular distribution of 
the particles. One can see that for reasonable values of /?rec ~ 0.2, 
amplification of factor of ~ 2 is expected. 



2.2 Acceleration time scale 

One acceleration cycle for a particle of energy 7 approximately 
equals the gyration time tg = Inymc^ I eBc at this energ}!^, where 
m is the mass of the particle and B the magnetic field strength of 
the upstream region. At the previous cycle the particle has A~^ less 
energy (and the cycle lasts proportionally shorter). The total time 
(acceleration time) for a particle of Lorentz factor 7 is the sum of 
the time spent in all the cycles until it reaches 7: 



^ The energy changes moderately during the cycle; we consider the final 
energy of the particle during the cycle slightly overestimating the gyration 
timescale and therefore the acceleration time. 
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Figure 2. Energy amplification A of particles for every Fermi acceleration 
cycle as function of the reconnection speed ySrec- The solid curve corre- 
sponds to particles crossing the reconnection layer at an angle 6 = and 
the dashed curve to an isotropic particle distribution. 
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For reasonable reconnection rates, the amplification is A ~ 2 which 
corresponds to acceleration taking place at a timescale of the order 
of the gyration period for the maximum energy of the particles. 

For a particle moving at a small angle with respect to the x di- 
rection {6 ~ njl case), the acceleration time scale can be estimated 
by setting the energy of the particle ymc^ equal to efirecBx, where 
X is the distance travelled along the layer. The acceleration time is, 
therefore, 4cc - ^/c = ymc^/(ej3recBc). For/?rec ~ 0.1 - 0.2, the 
previous expression results in face ~ ^g- This estimate is similar to 
that of eq. (5) showing that the 6 distribution of the particles has 
little eff'ect on the timescale for particle acceleration. 



3 CONSTRAINTS ON THE MAXIMUM PARTICLE 
ENERGY 

A number of physical processes can limit the maximum energy 
that a particle acquires. The acceleration ceases once the dynami- 
cal timescale of the system becomes comparable to the acceleration 
time 4cc or when the particle can no longer be confined within the 
reconnection region or when cooling of the particle is fast enough 
to inhibit further acceleration. Each of these processes is discussed 
in turn. 

In the following estimates, the magnetic field strength of the 
jet plays a critical role. It is determined by the (Poynting) luminos- 
ity L, bulk Lorentz factor F of the jet and the distance R from the 
central engine: 



B: 



Li/2 



(6) 



ci/2/?F' 

where the magnetic field B strength is measured in the rest frame 
of the jet. For observed energy £ of a UHECR, the energy of the 
particle in the rest frame of the jet is E' = E/T = ymc^. Through- 
out this work, we consider proton acceleration, i.e, m - m^ (the 
somewhat relaxed constraints for acceleration of heavier nuclei are 
straightforward to derive). 

Confinement Constraints: for the particle to be confined 
within the reconnection region, the gyroradius Rg cannot exceed 
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the region's thickness /redj Setting R^ - /rec and using eq. (|6]) I find 
that 






(7) 



During their acceleration drift along the x direction, the particles 
may eventually leave the reconnection region. For a given length of 
the reconnection layer 4 the maximum energy that can be achieved 
is Em = TejSrocBk = e/3rock(L/cy^^/R. For k ~ tec, the last expres- 
sion gives a maximal energy that is a factor y^j-ec smaller than E^^^. 
On the other hand, in the model that is discussed in the next sec- 
tion, /x is sufficiently larger than /rec and the escape along the layer 
is no more constraining for the maximum energy than the estimate 
ofeq. (7). 

Scattering processes can lead to particle motion perpendicular 
to the plane of Fig. 1 (along the reconnecting magnetic field; i.e., 
along the z axis). Escape of particles in the z direction will depend 
on the size of the layer 4 and the nature of the scattering process. 
Particle escape in this direction is likely a stochastic process and 
may (in connection to the amplification per cycle A) determine the 
distribution function of the accelerated particles but does not limit 
the maximum achieved energies. 

Cooling constraint: the accelerated protons cool through syn- 
chrotron emission and photo-pion production. Both of these pro- 
cesses can limit the energy of the UHECRs. In practise, however, 
it can be shown that cooling is dominated by the synchrotron emis- 
sion (Waxman 1995). The proton- synchrotron cooling timescale 

- 4^2^ can be set equal to the gyration time InElecYB to 



derive 



^syn 



Ine^B 



{m^c^f^ 



r=VS<'".'"''r"(^) 



-1/4 



(8) 



where eq. (6) is used in the last step of the derivation. 

Timescale constraint: The time available for the accelera- 
tion process to take place is limited by the dynamical (expansion) 
timescale of the jet. Equating the time R/Tc that takes for the jet to 
double its radius (in its rest frame) to the acceleration time (taking 
it equal to the gyration time; see Sect. 2.2), one finds 



77time _ 



IttY 



(9) 



In the next section, these constraints are applied to specific sources. 



3.1 GRB jets 

GRB jets have been proposed as sources of UHECRs by Milgrom 
& Usov (1995), Waxman (1995), and Vietri (1995) (see also Rieger 
& Duffy 2005; Murase et al. 2006; 2008). Using the observed long- 
duration GRB rate (e.g., Guetta, Piran & Waxman 2005) the energy 
release in y-rays from bursts at the local Universe can be estimated 
to be ~ lO'^'^erg/Mpc^yr. This is comparable to the rate of energy 
release required to power the observed UHECRs with E > 10^^ 
eV (Waxman 1995). GRBs are, thus, an energetically viable source 
provided that they can accelerate UHECRs with similar efficiency 
to which they produce y-rays. 

The (isotropic equivalent) y-ray energy of long-duration 



^ This constraint may be somewhat relaxed if the particles move predomi- 
nately along the reconnection layer {6 ~ nil) while been accelerated by the 
advective electric field. 
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GRBs is Ey ~ 10^^ erg while their duration is ~ 10s. This cor- 
responds to typical y-ray luminosity of Ly ~ 10^^ erg/s. During the 
active GRB phases, the luminosity of the jet is larger because of the 
likely moderate efficiency <50% in producing y-rays, and the qui- 
escent intervals in between ejection events. Here, I conservatively 
assume typical flow luminosity of L ~ 10^^ erg/s. Furthermore, 
compactness arguments (e.g. Lithwick & Sari 2001) bring the bulk 
Lorentz factor of the flow into the hundreds: F ~ 100 - 1000. 

Normahsing L = 10^2^52 erg/s, F = lO^^Fs.s and the dis- 
sipation radius (to be estimated below) io R - 10^^/?i3 cm, the 
cooling and timescale constraints JS]) and ([9]) read E^^ - 1 x 
^^^"^^l^lR^i^^'^i^ eV and E'^' = 9 X lO^^^f F-i eV respectively. 
The last expressions show that proton energies of £ ~ 10^^ eV 
are within reach if the magnetic reconnection takes place at large 
enough distance. 

The efficient dissipation of magnetic energy through recon- 
nection presumes that the jet contains regions with reversing mag- 
netic fields. These field reversals may be imprinted in the flow from 
the jet launching region (i.e. one deals with a non-axisymmetric ro- 
tator at the central engine) or are developed further out in the flow 
because of MHD instabilities. In the former case (i.e., that of an 
oblique rotator), the length scale of field reversals can be straight- 
forwardly estimated to be (in the lab frame) of order of the light- 
cylinder of the central engine 4^^ ^ cTq^, where Tce is the rota- 
tion period of the central engine. In the rest frame of the flow, the 
field reversal takes place on a scale /rec = F/J^^ = FcrcE (see, e.g., 
Drenkhahn 2002). In this picture, reconnection takes place over a 
range of distances but is completed at a distance R^^c for which the 
time scale of reconnection k^cl^c equals the expansion timescale 
RITc. Assuming a millisecond period rotator as the central engine 
of the GRB, the reconnection radius is (see Drenkhahn & Spruit 
2002 for a more detailed derivation) 



^rec — 



r^cTcE 



= 3 X 10' 



p2 J 

6-1 



(10) 



Using the last expression for the dissipation radius predicted 
by the reconnection model, one can write down the various con- 
straints (7), (8), and (9) as 



. X 10^' 



.6_iL 



1/2 



1^52 



L2.5 



eV, 



p5/2^1/2 

£'^" = 12x10^^^-^^^ 

^M - ^-^^ -^^ 1/4 1/2 

^52 ^-1 



eV, 



(11) 



E'^^ = 9 X 10 



rl/2 

19^52 



eV. 



Interestingly, for 6 ^ 0.15 the first and third constraints give the 
same result. Proton energies up to £ ~ 10^^ eV are plausible. 

In this model, UHECRs are accelerated at a distance R^^c (see 
eq. (10)) which is much shorter than the distance where the jet de- 
celerates interacting with the external medium. Adiabatic losses 
could, therefore, reduce the energy with which the particles es- 
cape. However, because of the efficient dissipation of Poynting 
ffux at /?rec, the magnetic field strength drops with distance steeper 
than R~^ for R > R^qc- It is, therefore, reasonable to assume that 
at least the most energetic particles become effectively decoupled 
from the jet and avoid adiabatic cooling. Note that in models for 
UHECR acceleration where the flow maintains a fixed fraction of 
Poynting flux, the magnetic field strength scales as 5 oc l/R (see 
eq. (6)). In those models, after the particle acceleration is completed 
Rg/ctdyn ^ (RB)'^ = const.; even very energetic particles remain 
coupled to the flow potentially suffering severe adiabatic losses. 



3.2 Jets in active galactic nuclei 

Powerful AGN jets can also potentially accelerate protons to UHE- 
CRs. Since we are not aware of any such jet with, say, L ~ 10^^ 
erg/s within the Greisen, Zatsepin Kuzmin (GZK) radius (Greisen 
1966; Zatsepin & Kuz'min 1966), one has to assume that such 
sources are transient (see Waxman & Loeb 2009 for discussion on 
observational constraints on the existence of sufficient number of 
AGN jets to power the observed UHECRs). A powerful AGN jet 
of L ~ 10"^^ erg/s and bulk F ~ 3 that contains field reversals on 
scales /rec ~ cTcE with Tqe ~ 10^ s (i.e. some 8 orders of magni- 
tude longer than that of the GRB central engine due to the larger 
mass of the compact object) can also satisfy the basic constraints 
for proton acceleration to the highest observed energies. For param- 
eters relevant for AGN jets, eq. (11) gives: 

6x\0''-^^ eV, 
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20 ^ 0.5 ^ 5 
^1/4 1/2 
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(12) 



9x10 



rl/2 

19^48 



eV. 



It remains, however, to be shown that such powerful jets are occa- 
sionally activated in nearby AGNs (see Farrar & Gruzinov 2009 for 
a proposed mechanism for such AGN ffares). 

Unlike the GRB flow that is expected to consist of protons 
(and maybe neutrons) because of the extreme temperatures and 
densities of the launching region, AGN jets may contain heavier 
nuclei such as iron. For iron composition of the UHECRs, the con- 
straints for acceleration on the luminosity of the source are signif- 
icantly relaxed (see, e.g., Pe'er, Murase, Meszaros 2009; Honda 
2009; Dermer & Rozaque 2010). Magnetic reconnection in less 
powerful AGN jets of L ~ 10^"^ - 10"^^ erg/s can accelerate iron 
up to energies of 10^^ eV. 



3.3 Injection energy 

The particles that enter the acceleration cycle discussed here, have 
passed through a preacceleration phase that makes their gyration 
radius Rg comparable to the thickness of the layer S. The required 
energy of the preacceleration phase depends on the details of the re- 
connection geometry but may be roughly estimated by the follow- 
ing considerations. Plasma of density pin enters the reconnection 
region with speed /?rec and leaves it with density pout at the Alfven 
speed /?A = VoTOTo^) through the layer: /?recPin/rec = TA/^APout^, 
where the length of the layer is assumed 4 ~ /rec- The compression 
Pout/Pin in the layer is model dependent. Assuming, for example, 
a Petschek-type geometry Pout/Pin ~ 2cr^/^ (Lyubarsky 2005), the 
thickness of the layer is (5 ~ lrec^/2cr. 

The downstream magnetic field B" (in the outflow frame) 
contributes only a fraction to the total pressure (the rest com- 
ing from hot particles). Pressure balance across the reconnection 
layer, therefore, constrains B'' < B (for thermally dominated down- 
stream B" «: B). In the downstream frame, the preacceleration 
energy E'^^^^^^ is given by setting 6 = Rg = E'^^^^^JeB" which 
yields ^pj-eacc ~ eB"6 < eB6 ~ eBlrec^/lcr. For typical GRB pa- 
rameters and o- = 100cr2, E'^^^^^^ < lO^'^e^^^f F-^o""^ eV. This 
is to be compared with the thermal energy in the downstream: 



^thermal ~ O-^m^C 



1/2^,^C2 



10^^cr2^^ eV. Substantial preacceleration 



above the thermal energy may therefore be required for the mecha- 
nism to operate. Tearing instabilities in the current layer are, how- 
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ever, expected to lead to fluctuations in the thickness of the layer 6 
and to regions where the preacceleration requirements are substan- 
tially reduced leading to the bulk of the particle injection. 



4 DISCUSSION 

Magnetic reconnection results in directly observed energetically 
substantial, non-thermal, high-energy particles during solar flares 
(Lin & Hudson 1971; Lin et al. 2003) or in Earth's magnetotail 
(Terasawa & Nishida 1976; Baker & Stone 1976; 0ieroset et al. 
2002). The puzzling particle acceleration measured by Voyager 1 
and Voyager 2 well beyond the termination shock of the solar wind 
(Stone et al. 2005; 2008) may take place at reconnection regions in 
the shocked solar wind (Lazarian & Opher 2009; Drake et al. 2010). 
Beyond our solar system, magnetic reconnection may operate in a 
variety of relativistic environments accelerating particles. Magnetic 
reconnection in the pulsar winds (Coroniti 1990; Lyubarsky & Kirk 
2001), the termination shock of pulsar winds (Lyubarsky 2003) or 
GRB flows (Thompson 2006) may result in particle acceleration 
and emission from these sources. If MHD jets contain field rever- 
sals on sufficiently small scale (of the order of the light cylinder), 
reconnection can efficiently power the GRB (Spruit et al. 2001; 
Drenkhahn & Spruit 2002; Giannios 2008) and AGN jet emission 
(Sikora et al. 2005; Giannios, Uzdensky & Begelman 2009, 2010; 
Nalewajko et al. 2010). 

Here, I present a model for UHECR acceleration in magnetic 
reconnection regions in strongly magnetized (AGN or GRB) jets. 
I argue that the converging magnetized plasma off'ers an ideal trap 
for particles to be accelerated to the highest observed energies. The 
mechanism discussed here assumes that reconnection is fast /?rec ~ 
(0.1 - 0.2)/?A but does not depend on, poorly known, reconnection 
physics. It is an efficient mechanism that accelerates particles on 
a timescale of order of their gyration period. Protons can reach ~ 
10^^ eV in GRB and luminous AGN jets. 
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